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The 0{a) improved Wilson quark action on the anisotropic lattice is investigated. We carry out numerical 
simulations in the quenched approximation at three values of lattice spacing (a^^ = 1-2 GeV) with the anisotropy 
^ = aa/a-r = 4, where aa and a-r are the spatial and the temporal lattice spacings, respectively. Using the 
dispersion relation of mesons, the bare anisotropy 7f in the quark action is numerically tuned below the charm 
quark mass region with the statistical accuracy of 1 % level. The systematic uncertainties in the calibration are 
examined and found to be under control in the continuum limit. Then we compute the light hadron masses and 
find that they are consistent with the result of the UKQCD Collaboration on the isotropic lattice. The effect of 
the uncertainty in the calibration on the hadron spectrum for physical quark masses is also found to be under 
control. 



1. Introduction 

The anisotropic lattice, which allows a finer 
lattice spacing in the temporal direction, is ex- 
pected to be useful for physics such as the spec- 
troscopy of exotic states, finite temperature QCD 
and the heavy quark physics. Since the mani- 
fest temporal-spatial axis interchange symmetry 
is absent, the anisotropy parameters of the action 
should be tuned by imposing the conditions with 
which the Lorentz invariance is satisfied for some 
physical observables. A lot of effort on anisotropic 
lattices has been devoted to charmonium systems 
|2|-^, but because of the computational cost for 
calibration, the light quark on the anisotropic has 
not been explored so far. 

In this report, we present our study on the 
0(a) improved Wilson quark action in the light 
quark region for the range of quark mass from 
the massless limit up to around the charm quark 
mass on the quenched anisotropic lattices with 
three lattice spacings, a^^ —1-2 GeV, at fixed 
renormalized anisotropy, ^ — a^/a-r — 4, where 
Ua and a-r are the spatial and temporal lattice 
spacings. First the bare anisotropy in the quark 
action is numerically tuned so that the renormal- 
ized fermionic anisotropy is equal to that of the 
gauge field by imposing the relativistic dispersion 



relation of mesons. The systematic uncertainties 
due to the anisotropy are examined. Then we 
compute the light hadron spectrum and examine 
how the uncertainty in the calibration affect the 
spectrum at the parameters of physical interest. 
More detailed discussions were presented in . 

2. Quark action and dispersion relation of 
free quark 

The quark action is the same as the Fermilab 
action |^ but defined on an anisotropic lattice as 
has been discussed in Rcf. In the hopping 
parameter form. 
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where k„ and the spatial and temporal 

hopping parameters, r is the Wilson parameter 
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and ce and cb are the clover coefficients. In prin- 
ciple for a given Kg-, the four parameters Ko-/ Kt, r, 
Ce and cb should be tuned so that Lorentz sym- 
metry holds up to discretization errors of 0{a?). 

In this work, we set the spatial Wilson param- 
eter as r = 1/^ and the clover coefficients as the 
tadpole-improved tree- level values, namely. 



CE = l/u„ul, CB^l/ul 



(3) 



and perform a nonperturbative calibration only 
for 7i?. The tadpole improvement [|| is achieved 
by rescaling the link variable as Ui{x) 
Ui{x)/ua and U^^x) — > Ui{x)/ur, with the mean- 
field values of the spatial and temporal link vari- 
ables, u„ and Ut, respectively. This is equiva- 
lent to redefining the hopping parameters with 
the tadpole-improved ones (with tilde) through 
Kg. = K^ju^ and Kt = k-r/ur. Wc define the 
anisotropy parameter as 7f = Kt /ka- 
For later convenience, we also introduce n 

- = ^-2(7;^ + 3r-4) =2(mo. + 4), (4) 

where moo- is the bare quark mass in spatial lat- 
tice units. 

The free quark propagator for the action 
satisfies the dispersion relation 



coshi?(p) = 1 



2(l + mo + i^p^) ' 



(5) 



where pi = — sin pi, pi — 2sin(pi/2), and 
mo = mQa-f^F is the bare quark mass in tem- 
poral lattice units. The rest mass Mi and the 
kinetic mass M2 are defined as Mi = £'(0) and 
I/M2 = £,'^cPE/dpl\p^o- One can tune the bare 
anisotropy parameter ^p so that the rest and ki- 
netic masses give the same values For small 
mo and with r = 1/^, this leads the expansion of 
in fTT-o as 



If' = [1 + l^nl] . 



(6) 



The mo dependence starts with the quadratic 
term for r = 1/^; therefore the dependence on the 
quark mass is small for sufficiently small mo- For 
example, let us consider the case of a"'^ =4 GeV, 
which corresponds to our coarsest lattice in the 
simulation. The charm quark mass corresponds 
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Figure 1. Dispersion relation of the free quark for 



to mo ~ 0.3 and at this value is different from 
^ by only 3%. Up to this quark mass region, one 
can expect that the difference of "fp from ^ will 
also be small in the numerical simulation. 

With our choice, r = 1/^, the action (||) leads 
to a smaller spatial Wilson term for a larger 
anisotropy f . The question is whether the con- 
tribution of the doubler eliminated by the Wil- 
son term becomes significant for practical value 
of ^. Figure |l] shows the dispersion relation 
for several values of toq in the case of ^ = 4. 
Towards the edge of Brillouin zone, the relative 
energy E{p) — E{0) rapidly decreases to ~ 0.1 
in temporal lattice units. On the other hand, 
typical energy scale of quarks inside hadrons are 
about Aqc d — 200 MeV for light and heavy-light 
hadrons. For our coarsest lattice, ~ 4.0 GeV 
leads to the doubler 's relative energy of about 
400 MeV. This seems large enough for eliminating 
the naive excited state contamination due to the 
doublers from the ground state signals, although 
more nontrivial doubler effect such as the shift of 
energy through the mixing of the ground states 
and the doubler states is not excluded. Our finest 
lattice has a^^ ~ 8 GeV which would be suffi- 
ciently large to avoid the doubler effect. In the 
case of heavy quarkonium, typical kinetic energy 
scale is mw^, where v is the quark velocity. For 
charmonium, v'^ ~ 0.3 and typical energy scale is 
about 300 MeV. Since this is not sufficiently less 
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than the doubler's relative energy on the lattice 
of a^^ ~ 4.0 GeV, finer lattices should be used 
for the charnionium system. 

3. Calibration procedures 

On an anisotropic lattice, one must tune the 
parameters so that the anisotropy of quark field, 
^F, equals that of the gauge field, £,g- 

^F{P,lG;'i,lF) = £.g{P,1g]h,1f) =C (7) 

Although and are in general functions of 
both of gauge parameters (/3, jg) and quark pa- 
rameters (k, 7f), on a quenched lattice one can 
determine Cg = C independently of k and jp, and 
then tune 7f so that a certain observable satis- 
fies the condition (Q). In this work, we define 
through the rclativistic dispersion relation of 
meson, 

E^{p)^m'+pyeF + 0[{p'n (8) 

for calibration. In the above expression, the en- 
ergy and mass E and m are in temporal lattice 
units while the momentum p is in spatial lattice 
units. converts the momentum into that in 
temporal lattice units. For finite lattice spacings, 
the above dispersion relation only holds up to the 
0[(p^)^] correction term. In the continuum limit, 
this higher order term in a would vanish and the 
rclativistic dispersion relation would be restored. 

We fit our numerical data for to the form 
Eq. (H) and obtain the value of £,f for each input 
value of bare anisotropy jf- Then we linearly 
interpolate i^^^ in terms of jf and find jp, the 
value of 7f for which £,f = £, holds. 

In order to estimate the systematic errors we 
also use the dispersion relation that corresponds 
to the lattice Klein-Gordon action 

cosh E{p) - cosh E{p = 0) = pV2C|:g- (9) 

The difference between these two calibration con- 
ditions shows the typical size of the lattice dis- 
cretization errors. Expanding this expression in 
a, ^KG is related to as 

Ug =^f[1 -my 12 + 0{a% (10) 

The same input jf gives a smaller value for S^kg 
than ^F, and therefore the tuned bare anisotropy 
"fp results in a larger value in the former case. 



Table 1 

Lattice parameters. The statistical uncertainty 
of Ut- is less than the last digit. 



/3 


7g 




size a^r ^ [GeV 






5.75 


3.072 


12'^ 


X 96 1.100( 6) 


0.7620(2) 


0.9871 


5.95 


3.1586 


16^ 


xl28 1.623( 9) 


0.7917(1) 


0.9891 


6.10 


3.2108 


20^ 


xl60 2.030(13) 


0.8059(1) 


0.9901 



4. Numerical Results of Calibration 

4.1. Simulation Parameters 

Numerical simulations are performed on three 
quenched lattices with the Wilson plaquette ac- 
tion at /? = 5.75, 5.95, and 6.10 and with the 
renormalized anisotropy ^ = 4. For the values 
of bare anisotropy of gauge field 7g, we adopt 
the relation of jg to ^ numerically determined 
by Klassen in one percent accuracy |^ . Table || 
summarizes the simulation parameters. The con- 
figurations are fixed to the Coulomb gauge, which 
is particularly useful for the smearing of hadron 
operators. 

The lattice cutoffs and the mean-field values of 
link variables in Table ^ are determined on the 
smaller lattices with half size in temporal extent 
for /3 — 5.75, 5.95, and otherwise with the same 
parameters, while at /3 = 6.10 the lattice size is 
16^ X 64. The lattice cutoffs are determined from 
the hadronic radius tq proposed in ||l0|| , by set- 
ting the physical value of vq as Tq ^ = 395 MeV. 
The mean- field values, Ua- and Ur, are obtained 
as the average of the link variables in the Lan- 
dau gauge, where the mean-field values are used 
self-consistently in the fixing condition ||^. 

4.2. Quark field calibration 

The calibration of the bare anisotropy jf in 
the quark action is performed as described in the 
last section. The tuned bare anisotropy parame- 
ter 7^,, at which ^f = ^ holds, is determined in 
the region from strange to charm quark masses 
using the dispersion relation of the pseudoscalar 
and vector mesons. The meson energies are ex- 
tracted from the meson correlators with momenta 
p — n{2n/ L), where L is the spatial lattice extent 
and n = (0,0,0), (1,0,0), (1,1,0), (1,1,1), and 
(2, 0, 0), where we take averages over rotationally 
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equivalent momenta. The energies are fitted to 
both linear and quadratic forms in in order 
to obtain in each channel, where the linear fit 
always uses three lowest momentum states. We 
adopt the result of quadratic fit at /3 = 6.f0, ex- 
cept for a few lightest quark cases in which the 
statistical fluctuation in large momentum states 
is severely large. At /? — 5.75 and 5.95, hnear 
fit is used in the whole quark mass region, since 
quadratic fits which include higher momentum 
states suffer from larger discretization errors. 

We measure S^p for two to four different values 
of 7i?. By linear interpolation we find the value 
of in each channel of PS and V mesons. We 
take the average of "f*p{PS) and ^*p{V) as the 
central value and use the difference to estimate 
one of the systematic errors. At each k, ^*p is 
obtained within f% statistical error, except for 
the lightest quark mass region. Figure |^ shows 
the K dependence of 'yp at (3 = 6.10. Similar 
behaviors are observed on other two lattices. To 
determine 7^ precisely in the light quark mass 
region and to extrapolate to the chiral limit, we 
fit ^p to the form 

^(m,) =Co + C2m2, (11) 
If 

where the quark mass in temporal units niq is de- 
fined as niq = — ^). The critical hopping 
parameter Kc are determined by a linear extrapo- 
lation of PS meson mass squared in terms of 1/k 
with two largest k's. The result of the fits is listed 
in Table ||, and represented by the solid line in 
Figure ^j. The linear fit in seems quite success- 
ful, and ^p at the chiral limit is close to the tree- 
level value, ^. The statistical uncertainty in 7^, is 
estimated to be of the order of 1% for the whole 
quark mass region. We also carry out a quadratic 
fit in TTig, as shown in Fig. ^ by the dashed line, 
and find that the difference in ')*p(mq — 0) from 
the linear fit in is at most 1%. This differ- 
ence should be considered as the systematic un- 
certainty in the chiral extrapolation. 

4.3. Uncertainties in calibration 

To examine the uncertainty in the calibration, 
the following analyses are also carried out. 

(i). We measure the difference between the ^*p^s 
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Figure 2. I/7I, vs 1/k at /3 = 6.10. The solid line 
shows the fit linear in while the dashed line 
represents the fit quadratic in rriq. 

Table 2 



Fit results for ^*p. 





Co 


C2 


hie 


5.75 


0.2558( 9) 


0.230(12) 


0.12640(5) 


5.95 


0.2490( 8) 


0.189(15) 


0.12592(8) 


6.10 


0.2479( 9) 


0.143(14) 


0.12558(9) 



for the pseudoscalar and vector mesons, which 
signals the 0{aa) systematic error. They tend to 
decrease towards smaller lattice spacing, and is 
already consistent with zero at /3 = 5.95. 

(ii) . To estimate the size of 0{a?) system- 
atic uncertainty, 7^ is also determined using 
the lattice Klein-Gordon type dispersion relation, 
Eq. (^. The results with two dispersion relations 
come closer to each other by decreasing the lat- 
tice spacing. The sizes of this uncertainty in the 
light quark mass region are 3%, 2%, and 1% for 
(3 — 5.75, 5.85, and 6.10, respectively. The be- 
havior in the large quark mass region is consistent 
with what is expected from Eq. ([l0|). 

(iii) . We obtain the response of meson masses 
to a change oi ^p at the fixed k. The effect of 
uncertainty in ^p on the meson masses is less than 
1%, if 7^1 is determined at this accuracy. This 
result is applicable to the relatively heavier quark 
mass region, such as < to^, and therefore in 
this region, the errors in the calibration are under 
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Table 3 

Quark parameters used in the hadron spec- 
troscopy. The numbers of configurations are 200 
at /3 = 5.75 and 100 at /? = 5.95 and 6.10. 



IF 



values of k 



5.75 3.909 0.1240, 0.1230, 0.1220, 0.1210 
5.95 4.016 0.1245, 0.1240, 0.1235, 0.1230 
6.10 4.034 0.1245, 0.1240, 0.1235, 0.1230 



control. The light quark mass region is examined 
in the next section. 

5. Light hadron spectroscopy 

5.1. Calculation of hadron spectrum 

Taking the central value of = 7fj ^^e light 
hadron masses are computed in the strange quark 
mass region rrig < mq < 2mg on the same lattices 
used in the calibration, while with smaller statis- 
tics. The parameters are listed in Table H. In this 
region, we regard that ruq is sufficiently small and 
adopt the value of in the massless limit. 

The quark propagators are smeared at the 
source with Gaussian smearing function with the 
deviation ~ 0.4 fm, in the Coulomb gauge. For 
baryons, two of the quarks are set to have degen- 
erate masses. 

In order to avoid the ambiguities in the defini- 
tion of the quark mass, we extrapolate the hadron 
masses to the chiral limit in terms of the pseu- 
doscalar meson mass squared, instead of 1/k. We 
assume the relation 

rn%g{mi,m2) ^ B{mi + m2); (12) 

so that — 2Bmi holds for the degenerate 

quark masses, mi — m2. Instead of (i=l,2), 
one can use mpg{mi,mi)'^ as the variable in the 
chiral extrapolation. For vector mesons and octet 
and decuplet baryons, we also use the linear rela- 
tions. The linear fit looks successful. 

5.2. Spectrum at physical quark masses 

We compare our hadron spectrum at the physi- 
cal quark masses with the result in the continuum 
limit by UKQCD Collaboration on an isotropic 
lattice [|ll[. We do not extrapolate our data to 
the continuum limit for lack of the number of dif- 
ferent /3's as well as the statistical accuracy. We 



adopt two different inputs to set the scale: one is 
hadronic radius and the other is the K* me- 
son mass, which were also adopted in |ll|] . From 
each of these scales the physical u, d and s quark 
masses are determined respectively. (We do not 
distinguish the u and d quark masses.) 

While we show only the result with the scale 
set by ro, similar result is obtained with the scale 
set by K* meson mass. The hadron masses ex- 
trapolated or interpolated to the physical points 
are shown in Figure |3[ together with the results 
of UKQCD in the continuum limit. In our data, 
differences between the results at f3 = 6.10 and 
5.95 are rather large compared with the difference 
between (3 = 5.95 and 5.75. This could be par- 
tially due to the different a dependence of 0{aa) 
and O(a^) lattice artifacts, and also due to the 
statistical fluctuation. Our results of the hadron 
masses seem to approach the continuum results 
by UKQCD Collaboration on an isotropic lattice. 
We also compare the parameter J, 



J 



mvdmv 



dm?pg 



(13) 



mv /•mps=mK* /mn 



which was introduced to probe the quenching ef- 
fect |l^. It is known that the quenched lattice 
simulation does not reproduce the experimental 
value, J = 0.48(2), and gives about 20% smaller 
value. Our results are consistent with those of 
UKQCD on isotropic lattices in the quenched ap- 
proximation, as shown in Fig. ^. 

5.3. Systematic errors of the spectrum 
from calibration 

To estimate the effect of the uncertainty of cal- 
ibration on the spectrum, we obtain the spec- 
trum at the same k's with slightly shifted bare 
anisotropy, 7^ = -\- 5^f- We set 5^p = 
0.1, which implies about 2.5 % shift of bare 
anisotropy. The difference between the masses 
with 7^ and ^*p is slightly amplified toward the 
chiral limit. Even for the lightest mass in each 
hadron species, the difference is at most around 
1%. This implies that the uncertainty of hadron 
masses at the physical {u,d) and s quark masses 
are about half the uncertainty m ^p. With the 
relativistic dispersion relation, ^p at = has 
been determined at each (3 within about 2% am- 



6 




0.0 0.5 1.0 



Figure 3. Hadron masses and the parameter J. 
The Oct is set by using vq. FiUed symbols are the 
results of UKQCD Collaboration on the isotropic 



lattices |11 



biguity: the statistical error of 1 % and the sys- 
tematic error of 1 % in the form of fit. Therefore 
there is 1 % level uncertainty in the hadron spec- 
trum due to the uncertainty in calibration. This 
feature makes the anisotropic lattice promising 
for future physical applications. 

6. Conclusion 

We studied the 0{a) improved quark action 
on the anisotropic lattice with anisotropy ^ — 
ac/a-r = 4. The bare anisotropy 7^,, with which 
^i? = ^ holds, is determined for the whole quark 
mass region below the charm quark mass, includ- 
ing the chiral limit, in 1 % statistical accuracy. In 
the massless limit, there is also about 1 % system- 
atic uncertainty in extrapolating 7^ to = 0. 
We estimate the typical sizes of 0{aa) and 0{a?) 
systematic uncertainties as to be 4% at /? = 5.75 



and smaller for larger (3. We then calculated the 
light hadron spectrum and found a consistent re- 
sult with previous work on an isotropic lattice by 
UKQCD. The relative errors in hadron spectrum 
for physical quark masses are half of those in 7^, . 
Since finite lattice spacing errors tend to vanish 
as a decreases, we expect to obtain the hadron 
spectrum in the continuum limit within 1% un- 
certainty due to the anisotropy. This encouraging 
results suggest that the anisotropic lattice would 
already be applicable to quantitative studies of a 
few percent accuracy. To achieve higher accuracy, 
nonperturbative tuning of the clover coefficients 
is required. 

The simulation was done on a NEC SX-5 at 
RCNP, Osaka University and a Hitachi SR8000 at 
KEK. H.M. thank JSPS for Young Scientists for 
financial support. T.O. is supported by Grants- 
in-Aid of the Japanese Ministry of Education (No. 
12640279). 
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